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The accretion of hot slowly rotating gas onto a supermassive black hole is considered. The 
important case where the velocities of turbulent pulsations at the Bondi radius rs are low, compared 
to the speed of sound Cs, is studied. Turbulence is probably responsible for the appearance of 
random average rotation. Although the angular momentum at rs is low, it gives rise to the 
centrifugal barrier at a depth rc — 1^/GMbh ^ rs, that hinders supersonic accretion. The 
numerical solution of the problem of hot gas accretion with finite angular momentum is found taking 
into account electron thermal conductivity and bremsstrahlung energy losses of two temperature 
plasma for density and temperature near Bondi radius similar to those observed in M87 galaxy. The 
saturation of the Spitzer thermal conductivity was also taken into account. The parameters of the 
saturated electron thermal conductivity were chosen similar to the parameters used in the numerical 
simulations of interaction of the strong laser beam radiation with plasma targets. These parameters 
are confirmed in the experiments. It is shown that joint action of electron thermal conductivity and 
free-free radiation leads to the effective cooling of accreting plasma and formation of the subsonic 
settling of accreting gas above the zone of a centrifugal barrier. A toroidal condensation and a hollow 
funnel that separates the torus from the black hole emerge near the barrier. The barrier divides the 
flow into two regions: (1) the settling zone with slow subKeplerian rotation and (2) the zone with 
rapid supersonic nearly Keplerian rotation. Existence of the centrifugal barrier leads to significant 
decrease of the accretion rate M in comparison with the critical Bondi solution for 7 = 5/3 for 
the same values of density and temperature of the hot gas near Bondi radius. Shear instabilities in 
the torus and related friction cause the gas to spread slowly along spirals in the equatorial plane in 
two directions. As a result, outer (r > rc) and inner (r < rc) disks are formed. The gas enters 
the immediate neighborhood of the black hole or the zone of the internal ADAF fiow along the 
accretion disk (r < rc). Since the angular momentum is conserved, the outer disk removes outward 
an excess of angular momentum along with part of the matter falling into the torus. It is possible, 
that such outer Keplerian disk was observed by Hubble Space Telescope around the nucleus of the 
M87 galaxy in the optical emission lines. We discuss shortly the characteristic times during which 
the accretion of the gas with developed turbulence should lead to the changes in the orientation of 
the torus, accretion disk and, possibly, of the jet. 
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INTRODUCTION 

The nearest giant elliptical galaxy M87 hosts a super- 
massive black hole (SMBH) with a mass of (3 — 6) x 
10^Ms„„ (Ford et al. 1994; Macchetto et al. 1997; Geb- 
hardt and Thomas 2009) and is filled with a hot gas that 
intensely radiates in X-rays. The ratio of the gravita- 
tional radius Vg = 2GMbh / ~ 10^^ cm to the distance 
of M87 from the Earth 17.9 Mpc (Gebhardt and Thomas 
2009) is known to correspond to an angular size of 4- 
8 ^as and to be only slightly inferior in this parameter 
to the black hole (BH) SGR A* at the center of our 
Galaxy. In contrast to SGR A*, the nucleus of M87 
is observable in the optical and ultraviolet bands and is 
actively investigated in all ranges of the electromagnetic 



spectrum. The X-ray luminosity Lx ^ 7 x 10 erg • s ^ 
(Di Matteo et al. 2003) of the central 100 pc in M87 is 
negligible compared to the Eddington one. At the same 
time, the jet in M87 is observed in the radio, optical, and 
X-ray bands, while the absence of a cooling flow in the 
hot gas of the central regions of M87 is indicative of its 
heating by the mechanical energy of the jet, shocks, and 
gas outflows maintained by BH activity. 

There are different opinions on the nature of this ac- 
tivity - in principle, it can be related to BH spin-down 
(Blandford and Znajek 1977), but the closeness (in or- 
der of magnitude) of the mechanical luminosity Lm ^ 
3 X 10''^ erg • s'^ (Churazov et al. 2002; Di Matteo et al. 
2003) and the power of the accretion energy release ex- 
pected in the picture of spherically symmetric Bondi ac- 
cretion assuming the accretion efficiency to be ~ 10% of 
Mbc^ (Owen et al. 2000; Young et al. 2002; Di Matteo 
et al. 2003) is suspicious. The point is that the Chan- 
dra X-ray telescope resolves the zone of the Bondi radius 
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rs ^ lO^^g, but the complex structure of the brightness 
distribution aUows the temperature and electron density 
of the fully ionized gas to be determined only on average 
over a zone of 15" — 20" whose size is more than an or- 
der of magnitude larger than the Bondi radius (Forman 
et al. 2007). The temperature and density along with the 
SMBH mass determine the Bondi accretion rate Mb- 

It is the assumption about such energy release during 
spherically symmetric Bondi accretion that is commonly 
made now in cosmological simulations of feedback from 
rapidly growing BHs in the nuclei of forming or merging 
galaxies on their growth and the outflow of matter in 
these galaxies (Springel ct al. 2005). This feedback is 
known to be very important for the evolution of galaxies. 

In this paper, we would like to recall that a low but 
finite angular momentum of the accreting gas in galaxies 
plays a very important role. Since the ratio rs/rg ^ 
(c/cs)^ ^ 10^ is very large, the presence of tangential gas 
velocities v^p ~ Cs{cs/c) '-^ 3 x 10~'^Cs ^ 2 km s~^ in the 
zone of the Bondi radius should lead to the appearance 
of a centrifugal barrier at rc > Vg and the stopping of 
a spherically symmetric flow (see, e.g., Kolykhalov and 
Sunyaev 1980, see also Limber, 1964; Fabian, Rees, 1995; 
Beskin and Malyshkin, 1996). In what follows, Cg is the 
speed of sound; Cs ^ 700 km s~^ at the Bondi radius. 

M87 has an anomalously low (even for elliptical galax- 
ies) " spin" that characterizes the rotation of the galaxys 
stellar component. SAURON data (Emsellem et al. 
2007) give only an upper limit for the spin, which, nev- 
ertheless, corresponds to rotation velocities much higher 
than 2 km s~^. The gas in the zone of the Bondi radius 
is unlikely to have fallen there as a result of the mass loss 
by stars. However, the hot gas near the Bondi radius 
is most likely subjected to turbulent motions that arise 
from the mergers of galaxies with M87 and the activity 
of its nucleus (for a discussion, see Appendix 1). It is 
these motions that impart a finite angular momentum 
to the gas that determines the flow pattern in the zone 

In the absence of mechanisms for efficient cooling of the 
accreting gas at r (10 — 100)rg, advection-dominated 
accretion flows (ADAFs) can emerge in principle below 
Tc- In addition, the angular momentum should some- 
how be withdrawn beyond rs upward - against the gas 
flow. ADAF-type solutions have been suggested long ago 
(Narayan and Yi 1994; Fabian and Rees 1995). They arc 
actively discussed as applied to both SGR A* and M87 
(Fabian and Rees 1995) and are distinguished by a low 
radiative cooling efficiency during accretion. Below, we 
will show that another regime of accretion works in ellip- 
tical galaxies similar to M87. 

The combined effects of "saturated" electron thermal 
conduction, deviation of the electron temperature from 
the ion one Tg <Ti, and bremsstrahlung near Tc are ca- 
pable of cooling the accreting plasma flow and contribut- 
ing to the establishment of delayed gas settling into the 



comparatively cold torus near the location of the cen- 
trifugal barrier. In the presence of magneto-rotational 
instability (Velikhov 1959; Balbus and Hawley 1991), the 
torus should spread out in both directions along the ra- 
dius and form a standard accretion disk that carries mat- 
ter to the BH and an outer accretion disk at r > rc re- 
sponsible for the removal of angular momentum into the 
zone beyond the Bondi radius (see Kolykhalov and Sun- 
yaev 1980). Numerical simulations show that the accre- 
tion rate for such a solution is only an order of magnitude 
lower than that in the critical regime of Bondi accretion 
for a gas with the adiabatic index 7 = 5/3 and depends 
on the initial angular momentum, i.e., on the position of 
the centrifugal barrier. As the angular momentum grows 
and Xc — Tc/tb changes from 0.01 to 0.04, the accretion 
rate in the considered mode increases from 8% to 17% of 
the accretion rate in the critical Bondi flow with 7 = 5/3 
(see the Figure at the end of paper). 

Our interest in this problem was heightened notice- 
ably by the assumption that the outer disk above the 
centrifugal barrier predicted by Kolykhalov and Sunyaev 
(1980) could be associated with the Keplerian disk ob- 
served by the Hubble Space Telescope by emission in the 
Ha, [O II], [N II], and other hues (Ford et al. 1994; 
Harms et al. 1994; Macchetto et al. 1997) at distances 
0.03 < r/vB < 1 from the BH. It is this disk that allowed 
the authors of the papers cited above to determine the 
BH mass in M87 just as was subsequently done for the 
BH in NGC 4258 from maser radio lines (Miyoshi et al. 
1995). 

The suggested picture of accretion has quite a few ob- 
servational consequences being discussed here and out- 
standing points on which we continue to work. First of 
all, the modiflcation of the solution by Kolykhalov and 
Sunyaev (1980) for the case of a disk that exists due to the 
removal of angular momentum and that is under the ex- 
ternal pressure of the atmosphere of hot gas settling from 
the zone near the Bondi radius to the centrifugal barrier 
is of considerable interest. Allowance for the heating of 
this disk outside by the electron thermal conduction of 
the hot gas is also important. It is this mechanism that is 
most likely responsible for the disk heating and that pro- 
vides its comparatively high luminosity in optical lines. 

The most important observational consequences to be 
discussed below are the following: 

(1) An enhanced brightness in X-rays and ultraviolet 
lines of the narrow zone of intense gas cooling above the 
torus near the centrifugal radius where ~ lO^'' — 10^'^ 
of the total mechanical luminosity of the M87 nucleus 
should be radiated away. This luminosity should not ex- 
ceed the X-ray luminosity of the central source, which 
(according to a possible overestimate) reaches Lx ^ 
7 X 10"° erg s'^ (Di Matteo et al. 2003). 

(2) The influence of turbulent motions in the zone with 
r > tb on the characteristic time in which the jet direc- 
tion changes (see Appendix 1). 
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(3) The emission from the disk above the centrifugal 
radius in ultraviolet lines through its heating by electron 
thermal conduction from the hot accreting atmosphere. 

(4) The effects from the eclipse of the central part of 
the counterjet by the disks being discussed. 

Unfortunately, the suggested picture cannot explain 
the high ratio of the mechanical luminosity of the M87 
nucleus to its bolometric luminosity in all ranges of the 
electromagnetic spectrum. At the very least, this re- 
quires assuming that the solution of a thin standard disk 
(Shakura and Sunyaev 1973) ceases to work in the zone 
with r < (20 - 100)r£, and transforms into ADAF. 

OVERALL PICTURE 
Spatial Structure 

The flow scheme is presented in Fig. 1. The Bondi 
sphere of radius rs bounds the zone of SMBH gravi- 
tational influence. Outside it, the temperature and den- 
sity asymptotically reach their values at infinity. Too and 
rioo. The subsonic atmosphere with the settling of mat- 
ter along the streamlines indicated by the dashed curves 
in Fig. 1 is located between the spheres with radii 
and Tc : 

= f/GAlBH = (Ma^lB)VB = XcTb, (1) 

where Ma<^|B is the Mach number relative to rotational 
velocity at radius vb- The term "atmosphere" is justi- 
fied, because, as we show here, the barrier reduces greatly 
the rate of accretion M. As a result, the latter is re- 
alized in the regime of slow settling with radial veloci- 
ties of the order of a few percent of the local speed of 
sound. Accordingly, the inertial effects and heat release 
through braking in radial velocity here are unimportant, 
while the rate M is appreciably lower than the maxi- 
mum possible one in a transonic (at the adiabatic index 
7 = 5/3), spherically symmetric Bondi flow. Note that 
such a transonic flow is often called critical. At 7 = 5/3, 
the critical flow reaches the speed of sound at point r = 
(Stanyukovich 1960; Kato ct al. 2008) and, hence, is 
transonic. 

The scale Vc (1) defines the location of the centrifugal 
barrier (CFB) and the centrifugal funnel; in what follows, 
I is the specific angular momentum. This scale is small 
{tc < rB) if the Mach number Ma^ls = (w¥'|b)/(cs|b) 
in rotation velocity at the Bondi radius is small. A typ- 
ical value of Maipls ~ 0.1 is determined by subsonic 
turbulence beyond the radius rs (sec Appendix 1). In 
the case of low Ma^^ls, the atmosphere occupies the ge- 
ometrically extended segment rc < r < rs- 

Above the CFB in Fig. 1, the dynamical influence 
from the presence of angular momentum in the settling 
flow is weak - the centrifugal force can then be neglected. 
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FIG. 1: Structure of the accretion flow in the presence of a low 
angular momentum. Rotation causes the flow (that is close 
to spherically symmetric settling above the centrifugal barrier 
r > rc) to be restructured qualitatively. Restructuring affects 
not only the regions near and below the barrier r < Tc but 
also above it r > r^- More speciflcally, a dense outer disk 
that threads the rarefied atmosphere and that transports the 
matter back, i.e., toward the settling flow, appears. 



Therefore, at r > the flow is radial with approxi- 
mately straight streamlines. At r ^ rapid rotation 
that accounts for a significant fraction of the Keplerian 
one begins to dominate. The flow is blocked by the CFB 
"wall" near which the streamlines curve greatly (see Fig. 
1). The minimum of the gravitational and centrifugal 
potentials, ^sum ~ ^GM /r + P /2r'^yi, forms the torus 
center in Fig. 1; in what follows, rcyi is the cylindrical 
radius equal to the distance to the polar axis. The sub- 
sonic accretion flow slowly sinks toward the minimum of 
the potential ^sum (see Appendix 2). 

The trajectories of individual particles are spirals run- 
ning over the surface that is formed by rotating the 
dashed curve (streamline) in Fig. 1 around the polar 
axis. The spiral pitch hsp in the meridional plane shown 
in Fig. 1 is equal to the distance traversed by a particle 
in one turn around the axis. It is defined by the velocity 
ratio v^l ^Jv^ ■\- Vg, where vg is the meridional velocity. 
The pitch hgp decreases rapidly as the minimum of ^sum 
is approached, because the rotation velocity begins to ex- 
ceed greatly the settling velocity. This is attributable to 
strong cooling and a sharp density rise toward the mini- 
mum. The ratio hgp/r (where r is the local radius) is 
great at r > rB, of the order of unity in the atmosphere 
at Ma;p|s ^0.1, and; as was said above, small near the 
barrier. 

A simplified scheme with a list of the structures main 
links, fiow separation at the torus center, and an hierar- 
chy of scales is presented in Figs. 2 and 3. For clarity. 
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FIG. 2: Scheme of physical processes: 1 - ultraviolet emission 
from the outer disk removing angular momentum, 2 - torus 
cooling by X-ray emission, 3 - optical thin-disk emission. The 
ADAF contribution to the emission is very small. 
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FIG. 3: Main scales: turbulence covers the outer segment 

from 2-5 kpc to rs; the Bondi radius is rs ~ lO^rg; 

Tc ^ 10^ rg is the barrier position; ta is the thin disk- 
ADAF transition; the radius of the last stable orbit around 
the SMBH is Zr^. 



Fig. 2 is not to scale: the region of the barrier is 
greatly magnified. Figure 3 shows the actual ratio of the 
sizes. In Fig. 3, the turbulent velocity field described 
in Appendix 1 occupies the region from an outer size 
of 2 — 5 kpc to the Bondi sphere ; the Bondi sphere 
tb ~ lO^Tg is located inside the turbulent field; the CFB 
Tc ~ Kf'Tg under the Bondi sphere is two orders of mag- 
nitude deeper (the typical condition Maj^js ~ 0.1 is as- 
sumed to be met); the thin disk - ADAF transition (point 
rA in Fig. 3) occurs somewhere within the centrifugal 
radius Tc] the radius of the last stable orbit around the 
SMBH is 3rc,; the jet base lies in a region of ~ lOr^. 

In Fig. 3, the X-ray luminosity of the hot gas from a 
volume equal to the volume of the Bondi sphere is ap- 



proximately three orders of magnitude lower than that 
from the barrier neighborhood. This, in turn, is three 
to four orders of magnitude lower than the mechanical 
luminosity. An energy of such an order of magnitude, 
first, is carried away by advection in the SMBH and, 
second, is transferred mechanically by the relativistic jet 
into the surrounding space. The thin disk - ADAF transi- 
tion (Narayan and Yi 1994) is caused by dissipativc heat 
release through viscous friction. The plasma does not 
cope with the emission of energy from this heat release 
due to a small density, because M ^ \Q~^MEdd under 
our conditions, where MecIcI is calculated from EEdd 
with a canonical coefficient of 0.1. In the ADAF regime, 
the fiow thickness above the equatorial plane is of the 
order of the radius, as is shown in Fig. 2, while the co- 
efficient of a -friction is great: a 1. The boundary of 
the disk - ADAF transition r = ta could be established 
if the energetically most weighty contribution from emis- 
sion 3 in Fig. 3, which, besides, differs spectrally from 
contributions 1 and 2, were detected and identified in 
observations. 

The centrifugal torus that "props up" the atmosphere 
from below is thick under weak cooling of the atmosphere 
and neighborhoods of the torus, with a "doughnut" ra- 
dius rtor of the order of Tc. As will be seen below, such 
is the situation in the absence of electron thermal con- 
duction, because the cooling time tsrs ~ 10^ yr through 
hot-plasma bremsstrahlung near tb is comparable to the 
characteristic times for the cooling flows in the central 
regions of galaxies and clusters. It is very long for the 
accretion problem under consideration, where the Keple- 
rian period IkIb at the Bondi radius is ~ 10^ yr. 

In the paper it is supposed that the stable stratification 
of the atmosphere in the SMBH gravity field prevents 
angular momentum transport through the atmosphere. 
Therefore this transport is going through the outer disk. 
The stable stratification is formed due to cooling of the 
bottom of the atmosphere near the torus. The cooling 
causes the entropy drop 



ds/dr > (2) 



with decrease of radius in the atmosphere. 

The detection of emission 1 in Fig. 2 allowed the pres- 
ence of a thin disk in M87 to be firmly established (Harms 
et al. 1994; Ford et al. 1994; Macchetto et al. 1997). In 
our model with an outer disk in Fig. 1, this emission is 
associated with the transition layer between a cold disk 
and a hot overlying atmosphere. The expenditure on 
emission is covered through transition-layer heating by 
an electron heat flow going from the atmosphere. 
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Torus Cooling and Angular Momentum Removal 
over the Region above Vc through a Barely Moving 
Atmosphere 

The picture changes quahtatively from the case with 
virial atmosphere in the presence of efficient torus cool- 
ing through electron thermal conduction and radiative 
heat transfer. The torus thickness rtor decreases greatly 
{rtor ^ re) and the doughnut becomes thin (see Fig. 
1). As has been said in the preceding section, the mat- 
ter sinks into the thin torus along the streamlines indi- 
cated by the dashed curves in Fig. 1 as it cools down. 
The matter in the torus is in rapid centrifugal rotation. 
Since the torus temperature is considerably lower than 
the virial temperature Tyir{rc), the rotation velocity 
Vipir^) ~ vxirc) exceeds significantly the speed of sound 
Cs{rc)- 

Under conditions of supersonic centrifugal rotation, 
the a friction becomes dynamically efficient. Due to 
friction, dense outer and inner disks of small thickness 
h <C r grow from the thin dense torus - the thin dough- 
nut is spread into a disk (see Fig. f). This is caused 
by strong mechanical coupling of differentially rotating 
rings with different radii inside the disk through friction. 
The disks are dense, because bremsstrahlung, which be- 
came efficient when the temperature dropped and the 
density rose exponentially in the torus, remains decisive 
in the disks as well. In the disks, the shear frequency 
dvip/dr = Vshear ~ and, hence, the friction stress 
Wrip is significant. As has been said above, for M87 the 
inner disk is restructured into ADAF at some distance 
TA within the radius due to a low cooling efficiency 
by radiation in this disk. 

In the torus, the downward-settling flow separates into 
two flows: one goes rightward into the outer disk, while 
the other goes leftward into the inner disk (sec Fig. 2). 
Thus, the ring of feeding of these disks is associated 
with the torus, as was assumed in the KolykhalovSun- 
yacv model (Kolykhalov and Sunyaev 1980). Through 
the atmosphere (toward the accretion flow!), the outer 
disk removes, first, the angular momentum of the matter 
settling through the atmosphere and, second, the angu- 
lar momentum of the inner-disk matter and the ADAF, 
while the inner disk, which is a standard one in the seg- 
ment rA<r<rc (Shakura and Sunyaev 1973), supplies 
the ADAF, the jet, and the SMBH with a mass flow. 



SYSTEM OF EQUATIONS 

Under conditions of the galaxy M87, the Knudsen 
number at the Bondi radius tb (estimated from r^) 
is small. Near and above r^, there arc enough collisions 
to equalize the electron, Te, and ion, T;, temperatures. 
However, for slow rotation at tb considered here, the 
CFB Tc (1) lies deep under the Bondi radius Tc <C r^. 



^ 1, X = r jvB^ Xc = Tc/rB- Accordingly, there 
is much heat being released as the matter sinks to this 
depth. The removal of such an amount of heat gives rise 
to a region of steep temperature gradients at a sufficient 
depth. At this depth, the Knudsen number increases, the 
collisions arc not enough, the thermal balance between 
electrons and ions is upset, and Tg becomes lower than 
Ti. Therefore, the system of our equations is written for 
a two-temperature plasma: 

dp/dt H- div(pz/) = 0, (3) 
pdv/dt + p{v-'^)v + Vp + pV<i> = Q, (4) 



p—^+p{v-V)ee+PeA\vv+(}L\-vq + e-a{Ti-Tc) =0, (5) 
at 



pde,/dt + p{v-V)ei + p.divw + a(T, - Te) = 0. (6) 

The system that consists of the continuity equation (3), 
the dynamical equation (4), and the thermal balances 
for electrons and ions (5) and (6) describes the flow from 
infinity to the CFB (accretion through the atmosphere). 
In the atmosphere, we neglect the a -viscosity in Eqs. 
(4) and (6). The molecular viscosity is low. The partial 
pressures and energies in Eqs. (5) and (6) are 

Pe = nfcsTe, Ce = {S /2)kBTe/ TUp, 



Pi = nkBT^, Ci = {3/2)kBTjmp. 

The total pressure and the gravitational potential in Eq. 
(4) are 

p = Pe+Pi, <^ = -GM/r. 

The electron heat flux q in balance (5) is significant. 
Let us represent it as a harmonic mean. 



Cehc 



1 



1 

Qlim 



|VT,| 



(7) 



This form allows the constraint on the flux q to be taken 
into account at high Knudsen numbers (see, e.g., Back 
ct al. 1996). Equation (7) is composed of the Spitzer 
expression 



1.84 X 10"^ 



In A 

and the limiting flux 



ri^/2[ei.gs-iK-icm-i] (8) 



kuTr 



/0.3^/% 



"pel 



(9) 



In (8) the temperature Te is in Kelvins, In A ~ 30 
40 is the Coulomb logarithm; in (9) we took Cg 
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^(10/3)fcBTe/TOp. The coefficient (ehc = 0.1-1 in (7) 
allows for the decrease in thermal conductivity due to 
the influence of the magnetic field associated with MHD 
turbulence (see, e.g., Zakamska and Narayan 2003). 

The constraint on the Spitzer thermal conductivity 
is taken into account using Eqs. (7), (8), and (9). It 
becomes necessary if the heat flux q calculated from 
Spitzers formula (Spitzer 1962) qsp = —nspT^ exceeds 
the vacuum flux nkBTey^KBT^Jrn^, where mg is the 
electron mass. This occurs in the region of steep temper- 
ature gradients when the characteristic length Tf./\T'J\ 
becomes comparable to the electron mean free path or 
the mean free path becomes of the order of the radius. 
The harmonic mean (7) is used to constrain the flux in 
the region of steep gradients. Due to the development of 
plasma instabilities, the limiting parameter / turns out 
to be small (Back et al. 1996). Wc will vary it within the 
commonly assumed range / = 0.05 — 0.3. At / = 0.05, 
the combination f 0.2,^/^ ^Jmp/rrif, in (9) is 0.35 and, ac- 
cordingly, qii,n = 0.35/r?c^. 

The radiative losses were calculated from the formulas 

^ — ^hrs ~t- G-compti 



Chrs ~ 2.1 X 10" 



-■compt 



AarksiTe - {Trad))£radne/meC, 



(10) 



(11) 



where eti^s and ecompt arc given in [erg s^^cm"'^], n 
and Ue are in [cm~'^] and ri = 0.1. Equation (5) is the 
energy balance between the change in electron internal 
energy in a Lagrangian particle, (i) the work for adiabatic 
expansion of matter pedivw under the electron pressure, 
(ii) the heat transfer by thermal conduction (7), (iii) the 
cooling through bremsstrahlung losses (10), and (iv) the 
Compton cooling (at Tg > 4(Trad)) due to the scattering 
of photons from the central X-ray source by electrons 
(11)-^ 

This central source is associated with the radiation 
from a region close to the BH. The role of electron-photon 
heat exchange in quasars and in Seyfcrt galactic nuclei is 
discussed in Levich and Sunyacv (1971). The formula for 
Compton energy exchange between electrons and pho- 
tons (11) includes the difference between the electron 
and photon temperatures; {Trad) ~ 1 keV is the effec- 
tive temperature of the photon gas, electrical neutrality 
Zrig = iii ~ n is assumed, the plasma is a fully ionized, 
hydrogen one Z = 1, ax = 6.652 x 10~^^ cm^ is the 
Thomson cross section, £rad is the radiation energy den- 
sity from the central source, L is the luminosity of this 
source. 

The expressions for electron-ion heat exchange are 



V = {'i/Z)V2^Znee^\nk/ml''^{kBTef'^. 

As has been said above, we assume that Z = 1. In (12) 
V is the electron-ion collision frequency (Spitzer 1962). 



SPHERICALLY SYMMETRIC APPROXIMATION 
OF A FLOW IN AN ATMOSPHERE 

This Section is devoted to derivation of non- 
dimensional steady-state equations from dimensional 
equations presented above. The system of non- 
dimensional steady-state equations is used for numerical 
integration. We omit this technical Section in the text 
for astro-ph. The full paper is published in Astronomy 
Letters, 2010, Vol. 36, No. 12, pp. 835-847. 



ACCRETION ATMOSPHERE IN THE GALAXY 

M87 

The boundary conditions that complement system 
(14)-(16) from previous Section are 

h{x = oo) = 1, ie(oo) = 1, U{oo) = 1, t'^{xc) = 0, (17) 

where normalized quantities n — n/rioo, te — 7e/7oo, 
ti — Ti/Too are used. The last of conditions (17) follows 
from the requirement that there be no electron heat ex- 
change between the torus and the disk (see Figs. 1 and 

2)- 

The unknown accretion rate M (13) enters into the 
fourth-order system (14)-(16) with four conditions (17) 
as a parameter. Therefore, the solution (n,te,ti) will 
depend on the non-dimensional parameter m : 



h{x]m), te{x;rri), ti{x;m), 



where 



r 

TB 



M 



Mrnaxil = 5/3)' 



1 . / 5 



Mb(7) = 47r[rB(7)]^PooCoo(7), 



rB 



GM 

Coo (7)^' 



Coo (7) 



' 27fcBTo, 



The set of such solutions forms a one-parameter family 
(the m runs through the family). From this family we 
will select the unique solution (the unique value of m) 
that corresponds to the maximum of the ratio 



a(Ti - Te) = S{me/mp)iynkB{Tt - Te), (12) 



(Tyir \c) / {Tq |c) 



(18) 
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of the virial temperature Tmr{r) at the CFB rc,Xc to 
the electron temperature in the atmosphere near the bar- 
rier. 

The minimum of the density ratio {pvir\c) / {p\c) (see 
Appendix 2) corresponds to the largest temperature ratio 
{Tvir\c) / iTe\c) (18). In an adiabatic atmosphere with- 
out any losses through electron thermal conduction (15) 
and radiative cooling (10), the virial density distribu- 
tion Pvir{i") is formed. The solution of system (14)-(17) 
that includes the thermal losses is given by the distribu- 
tion p(r). The torus is thinnest for the highest density 
(see Figs. 1 and 2). The inner and outer disks are then 
thinnest (among the solutions of the one-parameter fam- 
ily). At a large ratio of the densities of the outer disk 
and the atmosphere, the outer disk can exist for a long 
distances inside the hot rarefied atmosphere, despite the 
great difference in temperatures and rotation velocities. 
In this case, the dense outer disk threads the atmosphere 
and is retained up to the Bondi radius or even farther 
(Fig. 2). Thus, the angular momentum of the settling 
matter is removed outward. 

A typical example of integrating system (14)-(17) with 
the solution selection condition (18) is shown in Figs. 
4 and 5. This calculation was performed for rioo — 
0.15 cm-3. Too = IkeV, and M = 3.6 x lO^Msun- 
The integration was from point r — lOr^. The cool- 
ing through electron thermal conduction and radiative 
losses causes the atmosphere to be restructured qualita- 
tively compared to the adiabatic case. More specifically, 
in the lower part of the atmosphere adjacent to the CFB, 
the temperature decreases greatly, while the density rises 
exponentially. 

The dramatic increase in the "weight" of the atmo- 
sphere (the rise in p, see Fig. 5) maintained in the 
gravity field through an approximately radial pressure 
gradient is associated with cooling. The accumulation of 
great masses roughly follows a hydrostatic balance with 
a decreased plasma temperature. A very dramatic rise 
in density makes the bremsstrahlung losses (10) efficient. 
Thus, the matter heated by gravitational contraction in 
Fig. 4 cools down in two stages. Initially, a condensation 
is formed through electron thermal conduction and, sub- 
sequently, the mechanism of radiative losses is switched 
on in a high-density plasma, which increases the density 
even more. 

The temperature inversion shown in Fig. 4 takes place 
due to heavy bremsstrahlung losses in the lower part of 
the atmosphere. The electron heat flux q (7) changes its 
sign at the point of maximum of Te{r). Below the max- 
imum, the flux q heats (!) the plasma, partially com- 
pensating for the bremsstrahlung losses. The boundary 
condition T'^{xc) =0 (17) is met in a small neighborhood 
of point Xc, which is not seen in the scale of Fig. 4. The 
intermediate region between the Bondi radius and the 
lower part of the atmosphere is thermally most stressed 
in flux q. The flux q in this region is of the order of 




0.01 0.1 1 10 



FIG. 4: Distributions of the electron, Te, (solid curve) and 
ion, Ti, (short dashes) temperatures in the accreting atmo- 
sphere above the CFB at a typical value of Xc = 0.03 that 
follows from turbulent estimates in Appendix 1. The param- 
eter Xc (1) characterizes the rotation velocity at the Bondi 
radius. The values of T were normalized to the tempera- 
ture at infinity. For comparison, the long dashes indicate the 
variation in virial temperature Tt,ir. As we see, the virial 
atmosphere is much hotter in the lower part. 

the limiting flux (9), while the electron- ion collision fre- 
quency (12) is insufficient to equalize the temperatures 
Tg and Ti. 

In the intermediate region (narrow throat), the plasma 
is non- isothermal, Te < Ti (see Fig. 4). Above this re- 
gion, isothermality is maintained by an expansion of the 
area of the sphere 47rr^ through which the electron heat 
removal goes. In this case, the heat flux q per unit area 
of the sphere decreases. Below the throat, isothermality 
in Fig. 4 is restored by an exponential density rise. The 
point is that the limiting flux (9) is proportional to the 
density. 

The atmosphere is found to be essentially subsonic (see 
Fig. 5). The reduction in Mach number Ma at infinity is 
caused by an expansion of the sphere Anr'^ . In the lower 
part of the atmosphere. Ma decreases rapidly due to a 
sharp density rise and, accordingly, a rapid decrease in 
radial settling velocity according to integral (13). The 
maximum of Ma is associated with the intermediate re- 
gion. Since the radial velocities Vr are low compared to 
the speed of sound, the accretion rate M is also found 
to be low compared to the limiting one (13) which cor- 
responds to the critical Bondi flow with 7 = 5/3. The 
dependence of the dimensionless m (13) on the rotation 
velocity at the Bondi radius for typical rotation velocities 
is shown in Fig. 6. As we see, a significant reduction in 
accretion rate compared to the limiting possible one (13) 
is associated with the CFB. 

In conclusion, note that we considered the influence of 



8 



n / riip, 
10000 



1000 









U. 1 

Ma 


E 


\ / 


f \ 

\ 






\ / 
\ / 


\ 

^ Ma = V 


0708- 




\ / 


\ 




> 

n . 


\ / 

' \ 


\ 

ft 
\ 

\ 


0706- 


■■vir \ iv 


\ 


0704- 


.•A 

/ 


\ 

V \ 




i 1 
1 
1 

M 





0702- 

\ 

s 



0.01 0.1 1 x = r/rB 10 



FIG. 5: Sharp density rise in an atmosphere with coohng: cf. 
n(r) (sohd curve) for the solution of (14)-(18) with the virial 
distribution n„ir(r) (the curve with long dashes) . The arrow 
marks the enormous exponential density difference between 
n(r) and n^irir). Due to the presence of CFB, the matter 
settling is subsonic - see the dependence of the Mach number 
on radius (the curve with short dashes). Case Xc = 0.03. 



. dotm 































0.01 0.02 0.03 0.04 

FIG. 6: Influence of angular momentum on the accretion rate 
for typical (see Appendix 1) rotation velocities of the order of 
100 km s^^ at tb- The remaining parameters are the same 
as those in Figs. 4 and 5. 

electron thermal conduction on the structure of an ac- 
cretion flow with gas rotation in the case of a SMBH. 
We showed that the combined action of thermal conduc- 
tion and radiative cooling gives rise to a very massive 
thin cold torus and two disks, inner and outer. Due to 
heavy radiative losses, the inner disk near the torus and 
the outer disk all the way from the torus to the Bondi 
radius are found to be dense and thin. At smaller radii, 
the inner disk under the torus is restructured into a non- 
radiative ADAF. This is related to the conditions on the 



plasma density and temperature at the Bondi radius, un- 
der which a deficit of matter takes place. Because of this 
deficit, the emissivity of the inner disk at small radii un- 
der the torus is insufficient to maintain the regime of a 
standard disk by Shakura and Sunyaev (1973) (whence 
ADAF). The angular momentum of the accreting matter 
is removed over the outer disk. The outward transport of 
part of the matter beyond the Bondi radius is associated 
with the mechanical removal of angular momentum. This 
matter is in Kcplerian rotation. The possibility of out- 
ward transport of mass and angular momentum through 
the settling flow in the atmosphere is related to strong 
cooling and a high density of the matter being trans- 
ported. 

One of authors (R.A. Sunyaev) thanks Scott Trcmaine 
and Glenn van de Ven for discussing the problem concern- 
ing the "spin" of galaxy M87. Authors are very grateful 
to V. Astakhov for his careful translation. 
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